Spin gases as microscopic models for non-Markovian decoherence 



L. Hartmann 1 , J. Calsamiglia 1 , W. Diir 2 and H.-J. Briegel 1,2 
1 Institutftir Theoretische Physik, Universitat Innsbruck, Technikerstrafie 25, A-6020 Innsbruck, Austria 
Institutfur Quantenoptik und Quanteninformation der Osterreichischen Akademie der Wissenschaften, Innsbruck, Austria. 

(Dated: February 1, 2008) 

We analyze a microscopic decoherence model in which the total system is described as a spin gas. A spin gas 
consists of N classically moving particles with additional, interacting quantum degrees of freedom (e.g. spins). 
For various multipartite entangled probe states, we analyze the decoherence induced by interactions between the 
probe- and environmental spins in such spin gases. We can treat mesoscopic environments (~ 10 5 particles). 
We present results for a lattice gas, which could be realized by neutral atoms hopping in an optical lattice, and 
show the effects of non-Markovian and correlated noise, as well as finite size effects. 
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Decoherence is a fundamental mechanism believed to be 
responsible for the transition from the quantum to the classi- 
cal world Q|. Interactions between system degrees of freedom 
and (uncontrollable) environment degrees of freedom lead to 
entanglement, manifesting itself in the decoherence of the sys- 
tem state. Many decoherence models have been discussed 
in the literature, most prominent among them oscillator bath 
models |2| and spin bath models @|- Here, we study a phys- 
ically motivated model of a mesoscopic inhomogeneous spin 
bath. In particular, we describe the joint state of system and 
environment by a spin gas 0]. A spin gas is a system of 
quantum spins with stochastic time-dependent interactions. A 
physical model of a spin gas is a system of N classically mov- 
ing particles with additional, internal spin degrees of freedom. 
Upon collision, these quantum degrees of freedom interact ac- 
cording to some specified Hamiltonian. Hence, in such spin 
gases, classical kinematics drives the evolution of the quantum 
state, and also the decoherence of probe systems. In general, 
multiple non-consecutive collisions of particles are possible. 
In this sense, a spin gas provides a microscopic model for non- 
Markovian decoherence. 

In this letter, we determine an effective, time-dependent 
map that describes the decoherence process in a spin gas. For 
specific system-environment interactions, we can treat meso- 
scopic environments (w 10 5 particles) exactly and efficiently. 
To derive the map, we extend the description of certain states 
in terms of Valence Bond Solids (VBS) [ 5] to completely pos- 
itive maps. We do not restrict the analysis to single-qubit 
probes, but also consider the effect of decoherence on vari- 
ous (multipartite) entangled probe states (cf. jjj). Due to the 
stochastic nature of the interactions, our model does not dis- 
play symmetries, which otherwise simplify the treatment (see 
e.g. 0|] for homogeneous system-environment interactions). 
Throughout the paper, we concentrate on a specific realization 
of a spin gas, the spin lattice gas. However, our methods can 
be easily applied to decoherence in other spin gases, such as 
a Boltzmann gas |4|. Spin gases are not only toy models of 
theoretical interest, but could be experimentally realized with 
present-day technology and existing setups 0] . 

The model: The probe system A consists of Na qubits pre- 
pared in some arbitrary state. The qubits of system A interact 
with uncontrollable degrees of freedom of an environment B, 
leading to decoherence. We consider a microscopic decoher- 



ence model where the environment is described by a spin gas. 
The internal quantum degrees of freedom of this gas interact 
according to the time-dependent Hamiltonian 



H(t)=J29(r k (t),r l (t))H k i. 



(1) 



The function g depends on the physical nature of the pairwise 
spin interaction described by the Hamiltonian Hm- For such 
systems, we have shown Q |9[ that one can efficiently com- 
pute reduced density operators of up to ten particles even for 
mesoscopic system sizes (« 10 5 particles), if all Hamiltonians 
Hki commute and the initial state of the system is a pure prod- 
uct state. Here, we extend these results to take arbitrary initial 
system states and mixed environmental states into account, 
and thereby study the decoherence of multipartite probe states 
in a mesoscopic environment efficiently and exactly. 

For commuting Hamiltonians Hki, the joint state of system 
and environment at a time t, |* t ) = ]J k l e ~ lVkli - t ' )Hkl |* ), 
is determined by N(N — l)/2 interaction phases (fki(t) = 
Jo 9( x k{t'), xi(t'))dt' . These phases are associated with the 
adjacency matrix T(t) of a weighted graph G. The matrix el- 
ement Tki(t) = ifki(t) describes the neighborhood relation 
of particles k and I, or, equivalently, the interaction history. 
We are interested in the state of system A, i.e. the reduced 
density operator pA(t) = tr B\^t)(^t\- The commutation of 
the interaction Hamiltonians greatly simplifies the computa- 
tion of this operator: To compute this operator we first take 
into account only interactions between particles k E A and 
I E B, since interactions within B do not change the state of 
system A. In contrast to the general case of non-commuting 
Hamiltonians flfjll . entanglement in the environment does not 
influence the decoherence properties of the system. On top of 
that, the interactions within system A itself can be applied to 
the resulting state afterwards. 

In the following, we treat the case Hki — |H)ju(H| and 
initially completely polarized environment spins, \^)b = 
\+)® N where |±) oc |0) ± |1). Extensions to other commut- 
ing interaction Hamiltonians and arbitrary (possibly mixed) 
product environmental states are straightforward. The iso- 
morphism between completely positive maps and mixed states 
fTHl — together with a generalized Valence Bond Solids 
(VBS) picture Q U — determines the effective map &t 
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that takes an initial state at to = to the state at time 
t, i.e. PA(t) = £tPA(0)- This map £ t can be equiva- 
lently described by the state E t = 1 A ' <g> £f\<&) ($|. Here, 

1$) = ®l=i\0 + )k>k, oc |00) + |11), and A' is an 

auxiliary system with the same dimension as A. We can 
express the map in the Pauli basis a k , where <jq = 1, 
Then, £ t p = J2 ^k!...k NA ,h...i NA a k x ■ ■ ■ &k NA P^h ■ ■ ■ °i Na , 
where \ kl ...k NA ,h...i NA = (4>k 1 ...k NA \Et\4>i- L ...i NA ) with 

\4> kl ...k NA ) = 1 A ' <8> (cr fcl ...cr kNA ) A \$). The coefficients 
of the map are given by the coefficients of the state E t written 
in (tensor products of) Bell bases. As in the case of states 01, 
we can separately consider maps (or equivalently the states 
e[ 1 ') resulting from the interaction of the system with a single 
particle I in the environment. We find e[ 1) = 1/2(|$)($| + 

® k A li\Xk){Xk\) with \ Xk )k>k = 1/V2(|00) + e-^'Wlll)), 
where ip k i (t) is the effective interaction phase between par- 
ticles k £ A and I £ B. The state E t describing the total 
decoherence process incorporates the influence of all particles 
I £ B. We obtain E t (up to normalization) by calculating the 
Hadamard product of all E^ written in the standard basis, 
i.e. by component-wise multiplication. The matrix elements 
of Et expressed in the tensor Bell basis finally determine £ t . 
We find that E t has non-zero components only in the sub- 
space spanned by {\4 > k 1 ...k NA )} with kj £ {0,3}. The map 
£ thus contains only tensor products of Pauli operators 1 and 
a z . Equivalently, we can express the action of the map on an 
arbitrary probe input state p = J2 S s ' Pss' (0)|s) (s'| by deter- 
mining the evolution of the coherences p ss /(t). With s^jS^ 
we denote binary vectors of length Na- We can express the 
coherences as p BS '(t) = C ss i (t) p ss > (0) with |J] 

N B 

C ss >(t) = e^£*C^- s A)-r k JJ C0S [I( SA - s^) • T k ]. (2) 

fe=i 

The (rk)j = T k j for each particle k £ B are Na- 
dimensional vectors. The method described above can be con- 
sidered as an extension of the generalized VBS picture for 
states to one for completely positive maps. In this picture, 
we can determine the evolution of arbitrary system states in 
a mesoscopic spin environment, since the computational ef- 
fort to calculate the maps scales only linearly with the number 
of particles in the environment (as opposed to exponentially 
for general system-environment interactions). The size of the 
probe system is limited to about 10 spins for numerical com- 
putations due to exponential scaling with system size Na- 

The quantum properties of the system are directly linked 
to the classical statistical properties of the gas through T(t). 
In general, it is thus necessary to know the classical n-body 
phase-space distributions to give a complete description of the 
quantum state. If we assume no control over quantum or clas- 
sical degrees of freedom of the background gas we should av- 
erage over all possible collision patterns at any given time 0: 
C ail>l (t) = J dr>(r)C s v (r), where p t (T) is the probability 
that at time t the interaction history is given by V. 

For some gas models and regimes (like the Boltzmann gas 
studied in |4[]), correlations play a minor role and one can find 



analytical expressions for single-particle phase-space distribu- 
tions. In this paper, however, we study a lattice gas model that 
exhibits strong correlations, and produce the different random 
realizations of T(t) by direct simulation of the gas. 

The lattice model can be possibly implemented in a quan- 
tum optical system. It has already been demonstrated that an 
optical lattice can be used to store ultra-cold atomic gases. 
The degree of control in these experiments is extraordinary, 
opening the door to a wide ran ge o f ex peri ments and theoret- 
ical proposals ill W\ HI 111 l20l I2TII2I. One can choose 
a parameter regime where each lattice site is occupied by at 
most one atom The internal state of the atom (e.g. 

two meta-stable hyperfine states) can be stored in coherent 
superpositions over long time-scales (few minutes). Coher- 
ent inter-atomic interactions have been achieved by cold colli- 
sions fPUl . These correspond precisely to the Ising-type inter- 
actions chosen here. One can also find schemes |31 to induce 
a random (incoherent) hopping of atoms from one site to its 
neighboring sites. Hence, we consider anMxM lattice con- 
taining N particles that randomly hop from site to site with 
a hopping rate rj, and have nearest-neighbor interactions with 
coupling constant g D . With special relevance to possible ex- 
periments, we note that all results in this paper hold even when 
the environment particles themselves decohere. The only re- 
quirement is that the diagonal elements of the environment's 
state — in the canonical basis — remain unchanged. 

Decoherence of a single qubit: For a system consisting 
of a single qubit A = {1} and an arbitrary environment, the 
time dependent map corresponding to a particular collisional 
history is 

£tP = A 00 p + Xna- Z pa z + A01 {\po z - a z pl), (3) 

with Aoo = (1 + rcos7)/2, An = (1 — rcos7)/2, and 
A01 = O'r sin 7 )/2 where r(t) = n zeB cos(^), 7 (t) = 
J2ieB ^2~- Depending on the parameter regime, semi-quantal 
gases can follow various collision patterns. Accordingly, the 
dynamics of their quantum properties can differ considerably. 
If in every time step 5t a given particle collides with a dif- 
ferent particle and acquires an interaction phase 5 V , the dy- 
namics will be purely Markovian. The coherence of that par- 
ticle will decay exponentially fast with the number of time 
steps k = At/St, I poi| = [cos(^/2)] fc = e - A */ r = with 
r e Rj 8<5f/<5^. If, on the other hand, in a small time inter- 
val At a given gas particle has collided k times with the same 
particle, the coherent addition of the interaction phase leads to 
a Gaussian type of decay: |poi| = cos(kS v /2) sa e~ A * /( 2r s) 
with T g = 26t/S v . The exponential and the Gaussian de- 
cay are the two extreme cases, the dynamics of the coherence 
|poi I will usually lie in between. This also holds for the coher- 
ences in a multi-qubit density matrix. A complete character- 
ization of the system's decoherence is obtained by averaging 
the maps £ t over the possible collision patterns. The resulting 
map has the same form as l|3}, but replacing the coefficients 
Xij by their average values. An analytical expression for the 
time dependence of these averaged coefficients, and hence of 
the decoherence process, is in general hard to obtain. How- 
ever, for finite lattices, one can give a precise description of 
the time dependence. 
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Decoherence of bipartite entangled states: We now 

study the decoherence of different, initially entangled two- 
qubit states. In order to explore different regimes, we imagine 
a scenario where the probe particles can be displaced at a con- 
stant speed v relative to the gas. By varying the speed v and 
the distance d between the probe particles, we can highlight 
two effects: (i) By decreasing the probe speed, we analyze the 
effect of multiple interactions with the same particle in con- 
trast to interactions with different (independent) particles, (ii) 
By increasing the distance, we turn from correlated to inde- 
pendent collisions between probe and environment particles. 

(b) 




(a) 
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FIG. 1 : (a) Concurrence of two probe particles prepared in a max- 
imally entangled Bell state \<j> + ) oc |00) + |11> in a 100 x 3200 
lattice gas with N — 8 x 10 4 for g = O.877 and fixed probes 
(solid); and fast moving probes (dashed). The solid, light curve cor- 
responds to an analytical result obtained in a Markovian regime, (b) 
Concurrence of two probe particles at time t 1 = 25r]~~ 1 for differ- 
ent initial entangled states as a function of distance d between the 
probes : |^+) cx |01) + 1 10) (dashed-dark) , \(j>+) (solid-light), and 
|G)oc|+0) + (dotted). 

Figure^a) shows the decay of entanglement, measured by 
the concurrence Hal , of an initial Bell state \<p + ) in two ex- 
treme scenarios: (i) Fixed probe particles (v = 0). (ii) Large 
probe speeds, v/a ^> 77, where a is the inter-site spacing. A 
fixed value <p = 0.1 is assigned to the collisional phase every 
time a probe particle crosses an occupied site. These two sce- 
narios illustrate the difference between Markovian and non- 
Markovian environments. A large probe speed enforces a per- 
fect Markovian behavior which matches the analytical curve 
l24ll . Figure^b) shows the concurrence at a given time t\ as 
a function of the distance for three different entangled states: 
two Bell states and a cluster state (see figure caption). For 
Bell states the concurrence is equal to the absolute value of 
their only non-zero off-diagonal element in the density ma- 
trix, and therefore Figuref^provides direct information about 
the individual coherences. The figure shows the influence of 
correlated collisions: coherences poi,io ( m are robust 

against correlated noise, and coherences poo,ii ( m I0 + )) are 
especially fragile under correlated noise. The remaining co- 
herences decay in the same way under correlated or uncor- 
rected noise (hence, the weak distance dependence of \G)). 
From Fig. f^b) we also see that the immediate environments 
of each probe become more independent as d increases. 

The different behavior under correlated and uncorrected 
collisions can be readily understood. For two probe parti- 
cles (1 and 2) with very similar collision patterns, i.e., T%j = 
Tij+6j for all j, coherences associated with |01)(10| will only 
decay by a factor 2~ Nb J2 Sb e i5 ' SB , while |11)(00| will be 
"super-damped" by 2~ Nb J2 Sb e i{2T ^ +5 > SB . Therefore, clas- 
sical correlations in the collisions, e.g. induced by the geome- 
try of the set-up, can significantly influence the entanglement 
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properties of the system. 

Decoherence of multipartite entangled states: We now 

apply our method to investigate the decoherence of different 
multipartite entangled probe states of Na qubits. Due to the 
lack of simple, computable multipartite entanglement mea- 
sures, we use the negativity of bipartitions ffl as an indicator 
of multipartite entanglement in the system. That means, we 
consider bipartitions of the system, i.e. a partition consisting 
of a set of particles Ak and its complement Af., and investigate 
the entanglement properties with respect to the 2 Na ~ 1 — 1 in- 
dependent bipartitions tl2n . In general, we get a broad picture 
of multipartite entanglement in this way. For each bipartition 
we can determine its negativity A/A fc = f||o TAfc ||i — 1)/2 fl4ll . 
We define two multipartite entanglement measures: (i) the av- 
erage negativity Af, as the average over all bipartitions, AT 



!)E^^. and («)Mr 



i{A/aJ. 



Zero average negativity is a necessary condition for full sep- 
arability of the state, and N m i n = is a sufficient condition 
that the state is not multi-party distillably entangled. 




rjt r\t 
FIG. 2: (a) Log-linear plot showing the decay of the average negativ- 
ity M for six-qubit linear cluster state (solid), W-state (short-dashed), 
and GHZ states *) (dotted), |*') (dash-dotted), and |*") (long- 
dashed) in a lattice gas with 400 spins in a 40 x 40 lattice and cou- 
pling constant go = O.877. (b) Entanglement decay of four-qubit 
GHZ (long-dashed) and W (dotted) states in 20 x 20 lattice with 25 
spins and go = 8 x 10 -3 ?? illustrating finite size effects. Single- 
parameter (as) analytical fits are also shown (thin lines). 

We have examined different multipartite entangled probe 
states that interact with a lattice gas through a pairwise Hamil- 
tonian Hki = These states are (linear) cluster 

states |x) IH, GHZ-states and W-states \W) oc £^ \wi), 
where \wi) is the state corresponding to all spins in state 
zero except the i th that is in one. Since the noise process 
is basis dependent, we study variants of GHZ states corre- 
sponding to different local bases: |*) oc \0)® Na + \1)® Na , 

|*Vlo)|+) 0Ar - 4 - 1 + |i)h)® A '- 4 - 1 , |*")«l+)|o) 0Ar - 4_1 + 

| — ) \1)® Na ~ 1 . The decay of the average negativity is plotted 
in Fig. |3a). To qualitatively understand the different behav- 
ior of the curves we first derive analytic results in the limit 
of independent environments for each particle fc € A. In this 
limit, the decoherence process can be described by a tensor 
product of single-qubit dephasing maps S[ p = PkP + (1 — 
Pfe)4 p4* , where p k = ±(1 + l\ leB cos(2tp kl (t))). Co- 
herences decay as (t) = (2p k — l)poi(O). As before, 
the precise time dependence of decoherence will be given by 
the average p(t) = (pk)r(t) ov er different realizations. We 
assume that this average value is the same for every probe 



particle fc. Under the action of this dephasing map the famil 
of GHZ states remain diagonal in the GHZ-type basis 
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- 1 ~2)±\2 P -1\ Na ~ 1 } withA^ in = ±max{0,|2p-l| + 
l|^-i + |2p-l|^-l}. For VK-states, direct calcula- 



For such states, we can obtain the spectrum of the partial 
transposed operators with respect to any bipartition analyt- 
ically l23ll . and calculate average negativities for the three 
GHZ states given above. Here we give those with simple ex- 
pressions: Af = \\2p- 1\ Na , and A/ 7 ' = 2 n a- l -i Wn, 

\2p- 

tion leads to Af w = jj^w^ E^f 1 (^) y/a{N A -a). 

Several observations follow from these analytic results for 
independent environments: (i) Standard GHZ, |^), and W 
states remain A^-party distillable for all times, since M m in 
only reaches zero asymptotically as t — > oo. This does not 
hold for the | \&") GHZ, which has one partition that becomes 
disentangled at finite t, nor for |^'), for which Af vanishes at 
a finite time (all partitions are disentangled). In Fig. 2 we see 
that M also vanishes at finite time for the cluster state, (ii) In 
the limit of large system sizes, the average negativity is, to a 
very good approximation, given by the negativity of the half- 
half partition (the distribution of partitions with k particles is 
sharply peaked at k = Na/2). (iii) GHZ-type states decay 
exponentially as we increase the system size Na, while for 
W the coherences do not vary with the system size, rendering 
a weak dependence of its average negativity (constant to first 
order), (iv) For the states studied here, if a partition is initially 
more entangled than another, it will remain so also at later 
times. This does not hold in the presence of correlated colli- 
sions, which occur when the distance between probe particles 
is not large with respect to the relevant times t, i.e. d < ^/rji. 

In finite lattices and after long enough times the above de- 
scription fails due to inevitable correlated collisions. How- 
ever, the finite size leads to interesting effects that can also 
be easily understood. Figure 2(b) shows one such effect: the 
coherence (and hence the probe's entanglement) lost due to 
the stochastic interaction with environment spins, is partially 
recovered after a characteristic revival time. After long times 
a given environment particle will have collided nu times with 
a probe particle k. After s 3> M 2 steps, we describe the dis- 
tribution of values by a Gaussian of mean value (rife) = 



n = 4s/ Al 2 and the variance a 2 = ((rik — n) 2 ) = an, where 
a depends on the particular lattice model. We expect that at 
long times most blocking effects between environment parti- 
cles will be washed out, and therefore assume that different 
environment particles will have independent collision distri- 
butions. Hence, the total effect of Nb on a particular coher- 
ence will scale as C Nb , where C is the decay factor of the 
coherence due to a single environment particle. The value 
of C is given by an average taken over a Gaussian of mean 
<f and width a v : C = (cos(2</?)) = cos(2<,2 ) exp(— 8a 2 ). 
The phase tp is a sum of the collisional phases (with the cor- 
responding signs) involved in the particular coherence. For 
example, for a standard GHZ all phases are added tp = 
^2^=1 Pik leading to a mean value tp a m NatiS^ and a vari- 
ance a 2 ^ w a'NAnS 2 , l25ll . For the H^-state and a given co- 
herence, say p WlW2 , the phase is tp = <p\k — <P2k, leading to a 
vanishing mean value and to a purely exponential decay with 
a 2 ^ « a"2n<5^ — which is independent of the system size Na- 
The exponential decay of the VF-state and the periodic revival 
of the GHZ state can be clearly identified in Fig. 2(b). 

Summary: We have studied a microscopic, exact model 
for non-Markovian decoherence where the joint system is de- 
scribed by a spin gas. Using a generalized VBS picture for 
maps, we have determined the time-dependent maps for the 
decoherence process. We studied the decay of entanglement 
for different multipartite entangled probe states in a lattice gas 
with a mesoscopic number of particles. Depending on the pa- 
rameters of the gas, we have shown how to reach qualitatively 
different regimes such as Markovian- and non-Markovian, 
and correlated and non-correlated decoherence processes. For 
finite lattices we find that, although the interactions with the 
environment are stochastic, entanglement in the probes can 
spontaneously revive at a time given by the size of the lattice 
and independent of the number of environment particles. 
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